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Abstract 

First we describe briefly an information-action method for the 
study of stochastic dynamics of hamiltonian systems perturbed by 
thermal noise and chaotic instability. It is shown that, for the ensem- 
ble of possible paths between two configuration points, the action prin- 
ciple acquires a statistical form {5 A) = 0. The main objective of this 
paper is to prove that, via this information- action description, some 
quantum like uncertainty relations such as (A A) > -^=- for action, 

(Ax)(AP) > - for position and momentum, and (AH) (At) > -^=- 
for hamiltonian and time, can arise for stochastic dynamics of classi- 
cal hamiltonian systems. A corresponding commutation relation can 
also be found. These relations describe, through action or its conju- 
gate variables, the fluctuation of stochastic dynamics due to random 
perturbation characterized by the parameter rj. 

PACS numbers : 02.50.-r (Stochastic processes); 05. 40. -a (fluctuation); 
45.20.-d (classical mechanics); 46.15.Cc (Variational methods) 

1 Introduction 

The discussion of this work is limited to mechanical systems without con- 
sidering the quantum effect. We use the term regular dynamics to mean the 
mechanical processes with time reversible trajectories (geodesies) uniquely 
determined for each system by the Hamiltonian equations, or equivalently, 
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by least action principle. Compared to this deterministic character of regu- 
lar dynamics, one of the strong difference of irregular (random, stochastic, or 
statistical) dynamics is the uncertainty (unpredictability) of the trajectories 
of the system. This uncertainty can be illustrated in Figure 1 showing the 
diffusion of perfume molecules in the air. Clearly, from the point of view of 
classical mechanics, if there is no random perturbation from the molecules 
of air, an isolated perfume molecule leaving the hole a of the bottle at time 
t a has only one possible path or a fine bundle of paths having least action 
and arrives at a sole point b at time U. With the random perturbation of 
thermal noise of air molecules, however, the perfume molecules can arrive at 
many very different points b at a time t\>. This is the first dynamical uncer- 
tainty due to the random noise. In the theories of chaos, this uncertainty is 
usually measured by Kolmogorov- Sinai entropy or by Lyapunov exponent [1]. 
Another uncertainty is that, at time t&, the perfume molecules can arrive 
at a point b through different paths. This uncertainty in path space has 
been considered, e.g., in the path integral approach to quantum[2] and non 
quantum[3] dynamics and in the large deviation theory [4]. In general, these 
two uncertainty are not independent from each other. It is obvious that, 
without the first uncertainty, the second one cannot take place. 

While the mathematics of regular dynamics can be perfectly formulated 
on the basis of the least action principle in classical mechanics, irregular or 
random dynamics, found in diffusion, chaotic and other nonequilibrium phe- 
nomena, is much more complicated to be described due to its stochastic and 
random feature. Nowadays, statistical and thermodynamic theories of irreg- 
ular dynamics are still in development, among which we can cited nonequilib- 
rium statistical mechanics[5, 6], chaotic dynamics theory [1], anomalous trans- 
port theory [7], large deviation theory [1], small random perturbation theory [8] 
and path integral method starting from Brownian motion[3]. Many questions 
and topics concerning, e.g., time irreversibility [9], variational approaches[3, 
10, 11, 12, 13], chaotic nature of dynamic process[14, 15], connection with 
quantum physics[16, 17, 18] are still open to investigation. Much efforts have 
been made to clarify the origin of diffusion laws [19] such as the Fokker-Planck 
equation, the Fick's laws, the Fourier law, the Ohm's law and the anomalous 
diffusion laws (fractional or nonlinear) [7] . 

As a starting point of what we will describe in this paper, we would like to 
mention some characteristics shared by large deviation theory, perturbation 
theory, and path integral approach. In these theories, variational method 
is used to find the most probable (or optimal) paths (histories, trajectories) 
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with the help of rate functional [1], action functional[5, 8, 12, 13], or path 
integral[3, 18]. All these functionals can be called effective action S whose 
optimization, through a postulated exponential transition probability with 
a factor exp(— aS), allows one to find the most probable paths. Note that 
S are not necessarily the mechanistic action defined in classical mechanics 
with the Lagrangian L = E — U[3], where E and U are respectively kinetic 
and potential energy. Here L is defined for hamiltonian systems satisfying 
(in one dimensional space) 

• dH a p dH m 
X= dP andP = -^ (1) 

where x is the coordinates, P = mx the momenta, and H — E + U is 
the hamiltonian of the system. The mechanistic action is defined by A = 
f% L(x, P, t)dt. Its stationary 5 A = according to least action principle leads 
to the Euler-Lagrange equations [20] 

d_dL_9L_ Q (2) 

dt dx dx 

If we consider the Legendre transformation H = Px—L, Eqs.(l) can be easily 
derived from Eq.(2). It is worth noticing that, when the effective actions 
mentioned above are calculated by time integral of an effective Lagrangian, 
Eq.(2) is always satisfied by the optimal paths[3, 5, 8, 13, 18]. 

A point to be noticed is that, in the above theories and their applications, 
one is most interested by the most probable paths whose neighborhood pro- 
vides the basic contribution to the transition probability [3, 8, 12, 18]. The 
other less probable paths of larger actions are often neglected or only taken 
into account (or buried) in the path integrals[2, 3]. As a matter of fact, the 
application of classical action principle only to the optimal paths is, to our 
opinion, incomplete. The physics represented by the larger action paths is an 
inseparable part of the dynamics and may be essential for the fundamental 
understanding of irregular dynamics. Recently, an informational method was 
proposed [21, 22] to treat all the possible paths as an ensemble. The method 
considers hamiltonian systems under random perturbation of thermal noise 
and chaotic instability leading to the uncertainties shown in Fig. 1. These 
uncertainties can be measured by a path information associated with differ- 
ent paths between two points in configuration space. The path information 
is optimized in connection with the average action of the hamiltonian system 
calculated over all the possible paths. We obtain an exponential probability 



3 



distribution of action. In principle, this action can be any effective action 
mentioned above. But in our previous work, the classical mechanistic action 
is used. It is worth noticing that this variational method is proven to be 
equivalent to an "extended least action principle" for stochastic dynamics, 
i.e., instead of SA = 0, we have (SA) = where the average (•) is taken over 
all the paths. Using classical action turns out to be an useful choice since 
we have been able to derive the diffusion laws mentioned above for all the 
possible paths not only for the optimal paths. In this paper we will describe 
some mathematical consequences of this approach. Our main objective here 
is to derive non quantum uncertainty relations, which are in fact a necessary 
consequence of the description of stochastic dynamics with path probability 
given by distribution of action. 



2 Optimizing information- act ion 

First let us recapitulate briefly the method. Let pk(b\a) be the transition 
probability along a path k (k — 1,2, ...w) from a position a and to another 
position b. The dynamic uncertainty associated with pk(b\a) is measured 
with the following Shannon information 

w 

lab = ~^2Pk(b\a) \wp k (b\a). (3) 

fc=i 

We have the following normalization J2k=i Pk{b\a) = 1. If if is very large in an 
ergodic phase space and if the paths are sufficiently smooth, J2k should be re- 
placed by a path integral in the Feynman sense [2], e.g., I a b = — J D(x)pk(b\a) \np k (b\a) 
in keeping a and b fixed. The average action between a and b is given by 

w 

(A ab ) = Y,Pk(b\a)A ab (k). (4) 
k=i 

where A ab (k) is the classical action along a path k. Our optimal information- 
action method consists in the following operation: 

It! V) 

S[Iab + aY,Pk(b\a) - r)Y,Pk(b\a)A ab (k)] = (5) 

k=l k=l 

leading to 

Pk(b\a) = - exp[-r}A ab (k)], (6) 
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where the partition function Z = J2k exp[— rjA ab (k)} = J V{x) exp[—r]A ab (k)]. 
It is straightforward to see the following relationships : 



I ab = In Z + r](A ab ) 



(7) 



and 



(A ab ) 




(8) 



It is proved that [21] the distribution Eq.(6) is stable with respect to the 
fluctuation of action. If one uses the action of a free particle[2, 3, 6], pk(b\a) 
is just the transition probability of Brownian particles. In this case, we have 
a precise physical meaning of the multiplier r], i.e., T) = where m is the 
mass and D the diffusion constant of the Brownian particle. The signification 
of T) will be discussed later in a general way for particles moving in a potential 
field U(x). 

3 Extended action principle for irregular dy- 
namics 

As expected, Eq.(6) is a least action distribution, i.e., the most probable 
paths are just the paths of least action 5A ab (k) = satisfying Euler-Lagrange 
equation and Hamiltonian equations. The other paths do not satisfy Eqs.(l) 
and (2). In general, the paths have neither 5A ab {k) = nor 5(A ab ) = 0. 
Eq.(5) implies following relationship 



By using Eqs.(3), (4) and the distribution (6), it is easy to calculate that 
Eq.(9) is equivalent to 



On the other hand, in mimicking equilibrium thermodynamics, we can 
define a dynamic potential 



- rj5A ab + 6I ab = 0. 



(9) 



(5A ab (k)) = Y,Pk(b\a)6A ab (k) = 0. 



(10) 



k 



-InZ 
V 
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as an analog of the free energy of Helmholtz. From Eqs.(7) and (9), it is 
straightforward to see that the extended action principle Eq.(lO) is equivalent 
to the following variational principle 

5* = 0. (11) 

A remarkable application of this above variational approach to computation 
of thermodynamic properties of hamiltonian systems was (independently) 
carried out recently [10]. The authors derived thermodynamic equations of 
motions for equilibrium systems with different Hamiltonians, already known 
in the literature and used in simulations of molecular dynamics, by consid- 
ering the particle histories in phase space on constant energy surface. 

The extension of action principle has some consequences on the equations 
of motion as discussed in [23]. For example, for the paths whose action is 
not at stationary, we get &L q^ — 9L q^ an d P ^ — fjf which implies a 
stochastic equation like 

R (12) 



OX 

where R is a random force representing the random perturbation of thermal 
noise and chaotic instability. On the other hand, using the action principle 
Eq.(10), we recover Eqs.(l) and (2) for the ensemble of paths: 

lddL k (t)\ ldL k (t) 



dt dx \ dx 



(13) 



and 

and (P) = - ( ^ ) . (14) 





This implies that the mean of the "random force" R over all possible paths 
must vanish, i.e., (R) = 0, required by the extended action principle. The 
second equality of Eqs.(14) is the random analog of Newton's law and has 
the same content as the Feynman-Hibbs quantum Newton's law in Eq.(7-42) 
of [2]. 

In what follows, we are concerned with the spreads of the distribution 
of the stochastic dynamics, in other words, the deviation of the irregular 
dynamics from the regular one due to random perturbation. This uncertainty 
is a priori measured by the path information we introduced. But here it will 
be analyzed at the level of mechanical quantity like position, momentum, 
action and energy. 



6 



4 Uncertainty relations 



For the sake of simplicity, we suppose the point b at x is very close to the 
point a at xq and the transition takes place in the infinitesimal time interval 
St — t — to- This segment can be considered as a factor (a small element of 
a long path) in path integral technique[2]. Let 5x = x{t) — xo(t ), the action 
is given by 



A ab (x) = ^£ + F^5t - U(x )6t, 



(15) 



where F 



- ( lr ) is the force on the path and m the mass of the 

\dxj ( x+XQ )/2 v 

studied system. The transition probability is 



Pk(b\a) 



exp —7] 



rn 

m 



6x 2 



with 



exp 




771 r 2 c- 

2Si 5x+F J 6x 



(16) 



(17) 



2n5t 
mr) 



The potential energy of the point xq disappears after normalization because 
it does not depend on x. F is considered constant on small Sx. It is easy to 
show[22] that the probability of Eq.(16) satisfies the Fokker-Planck equation 
and that other diffusion laws can be trivially derived from it. 

How far are the randomly perturbed paths deviated from the optimal 
paths of regular dynamics? How different are Eqs.(14) from Eqs.(l)? This 
question can be answered, under different angles, by the standard deviation 
of action (AA) 2 = (A 2 ) - (A) 2 = (A 2 ) - (A ab ) 2 = -^f 1 , of position 
(Ax) 2 = (x 2 ) — (x) 2 , of momentum (AP) 2 = (P 2 ) — (P) 2 and of Hamiltonian 
(AH) 2 = (H 2 ) — (H) 2 . Using the formula below 1 , we obtain 



(X - Xq) 



F5t 2 

2m 



(18) 



1 J^° oo dxexp(—ax 2 + 2 r yx) 
i ^^[ 7 exp( 7 2 /a)]. 



exp(7 2 /a)y^ and x n dx exp(— ax 2 + 2^x) 



2 n ~ 1 a. V a c/7 71 
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(19) 



and 



N2v St F 2 5t 4 , , 

(( x - Xo f) = — + j 20 



Since (Ax) = (A(x — xq)), the above relationships imply 



(Ax) 2 = — (22) 
mrj 



and 



<AP>> = ^. (23) 

Cancelling the time St and m in the above two equations, we get a "classical 
uncertainty relation" for irregular dynamics 

(Ax) 2 (AP) 2 = \. (24) 

Notice that this is only an asymptotic relation for St — > 0. For a measurable 
(longer) length of time and path, we divide them into small segments of the 
order of St and 5x. Eq.(24) should be valid for each segment. Thus it can be 
trivially proven that the total deviation (A) 2 of whatever quantity on a path 
is a sum of all the deviations on its small segments (Gaussian law of errors, 
see for example [6]). So in general, we must write 

(Ax) 2 (AP) 2 > 1 (25) 

which is valid for any measurable period of time and path length. If the 
system is in rotation with Sx/R = 58 where 9 is the rotation angle and 1/R 
the curvature of Sx, then we get 

(Afl) 2 (AJ) 2 > 1 (26) 



s 



where J is the angular momentum. 

The action uncertainty can be calculated from Eqs.(8) and (17) : 

. . 4 1 F 2 St 3 

^> = 2S-^r (27) 



so that 

(AAY< 



d(A ab ) 
dr) 

1 

2^ 



(28) 



This relation can be written as (AL) 2 (St) 2 = where (AL) 2 = (L 2 ) — (L) 2 
is the standard deviation of Lagrangian. For arbitrary time and path length, 
we must write 

(AA) 2 > JL. (29) 

and 

(AL) 2 (At) 2 > JL (30) 

where (At) 2 = (t 2 ) — (t) 2 is the standard deviation of time measure. Eq.(30) 
can be verified in the same way as for position and momentum, through the 
calculation of (t) and (t 2 ) with the distribution Eq.(16) in relaxing St and 
fixing Sx. 

Now let H(t) be the Hamiltonian on the considered path segment at mo- 
ment t, we have the Legendre transformation H{t) = Px — L(t). Considering 
L = E — U and E = Pi/2, we get 

(AH) 2 = (AL) 2 + 4({EU) - (E)(U)). (31) 

E and U are two independent variable, so (EU) = (E)(U), we obtain 

(AH) 2 = (AL) 2 (32) 

and 

(AH) 2 (At) 2 > JL (33) 
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Before giving an interpretation of these non quantum uncertainty re- 
lations of irregular dynamics, i.e., Eq.(25), Eq.(26), Eq.(29), Eq.(30) and 
Eq.(33), we would like to indicate two points. First, it seems that Eq.(29), 
the action uncertainty relation, is the most essential one because it tell us 
the spread of the action distribution and the deviation of the extended action 
principle (SA) = from the conventional one 5A = 0. The other relations 
concerns the conjugate variables of action and can be a priori derived from ac- 
tion uncertainty. Second, the coefficient rj plays a central role in this descrip- 
tion of the dynamics. Roughly speaking, 77 — » represents large deviation 
of the perturbed dynamics from the regular one (large uncertainty), large 77 
represents small random perturbation and deviation (small uncertainty), and 
77 — > 00 implies vanishing perturbation and convergence of irregular dynamics 
to regular one (zero uncertainty). 

A quantitative relationship between the uncertainty and thermal noise 
can be shown with a special example, Brownian motion (F = 0) with diffu- 
sion constant D = /z/cbT[6], where fi is the mobility of the particles, ks is 
the Boltzmann constant and T the temperature. Here we see a relationship 
between 77 and temperature characterizing the thermal noise: 77 = 2m ilk B T ' 
In this case, we can write, for example, 



which implies that the paths may be distributed very widely in phase space 
around the optimal (least action) paths at high temperature. 

5 A commutation relation 

Now let us define a certain functional Gk{x) of paths. The average of the 
functional over all the paths is given by 



Now differentiate this equation with respect to x, the average does not depend 
on position, so d^W) _ q ^he derivative of right hand side yields two 



terms: EfcPfc(%) £ and -r]J2kPk(b\a)G k (x) °^ ' . This implies 



(Ax) (AP) > 2mfik B T, 



(34) 




(35) 




(36) 
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which is the classical analog of the Feynman-Hibbs relation given by Eq.(7- 
30) of [2]. Imposing Gkix) = x, we obtain[2]: 

\ m —&i— xi ) - \ x * m ^^/ = V (37) 

where % is the index of a segment for the time period St = ti — of a certain 
path k. This relation can be written as 

(P i+1 Xi) - (xiPi) = ~, (38) 

which can be considered as a classical commutation relation, where Xj is the 
position at time U, P,, + i = Xi+1 s ^ Xi and Pj = Xi ~^~ 1 are the momenta at time 
t i+1 and ti, respectively 

We would like to indicate that the above results, i.e., the uncertainty and 
commutation relations, were obtained by using the mechanical Lagrangian 
and action. Nevertheless, as mentioned above, the same conclusions can be 
reached with other action functionals and Lagrangians. It can be trivially 
shown that similar uncertainty relations arise with, e.g., Onsager-Machlup 
Lagrangian L OM = (x — /iF) 2 [3, 5] and Freidlin-Wentzell Lagrangian L m = 
[ip — b({p)} 2 where (p = b(<p) is the differential equation of some continuous 
function if of paths when the perturbation vanishes[8]. If the Lagrangians 
do not have mechanistic interpretation, the uncertainty relations should be 
understood differently depending on the physical content of the variables. 

6 Concluding remarks 

We have briefly presented an information-action method for stochastic dy- 
namics of hamiltonian systems. We see that in this approach the Euler- 
Lagrange equations, the Hamiltonian equations and the action principle only 
hold in averaged form over the ensemble of all the possible paths between two 
state points. The main result of this paper is to show that some uncertainty 
relations, very similar to those in quantum mechanics, exist for stochastic 
dynamics of hamiltonian systems. 

We are in the realm of classical physics, no quantum effect is considered. 
In addition, the systems under consideration are possibly macroscopic. There 
is normally the consensus that exact value of position and speed can be a 
priori assigned to a particle at the same time. Then how to understand 
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the uncertainty relations in classical dynamics? A plausible understanding 
is that (AA) > implies the distribution of action cannot be arbitrarily 
narrow around the least action, or the fluctuation of action has a non zero 
minimum due to the random perturbation. In the same way, (Ax) (AP) > - 
means that the fluctuations or the distributions of position and momentum 
cannot be simultaneously arbitrarily narrow. It has no meaning about the 
precision of the measure, because a classical body does have exact position 
and speed and the means (x) and (P) can be made arbitrarily close to the real 
values of the body whatever the fluctuation (Ax) and (AP) of the separately 
measured values. 

A question naturally arises. Do these classical uncertainty relations have 
something to do with their quantum analogs? Remember that in this work we 
consider classical systems under the random perturbation of thermal noise. 
This essential point is illustrated by the example of Brownian motion in 
Eq.(34) which tell us that if the thermal noise vanishes for T — > 0, we get, 
e.g., (Ax)(AP) > 0. So no quantum effect can be observed here. One 
needs additional hypothesis in order to enter into the quantal world from 
the present results. For example, if one assumes n = 1/h, the distribu- 
tion function Eq.(6) becomes an universal Brownian distribution of Fenyes- 
Nelson[16, 24] satisfying Fokker-Planck equation[23] whose forward and back- 
ward versions can lead to Schrdinger equation[16, 18]. At the same time the 
classical uncertainty relations become the Heisenberg's relations. To our 
opinion, the passage from classical stochastic dynamics to the stochastic 
(quantum) mechanics of Fenyes-Nelson is not automatic. The assumptions 
of an universal noise and the emergence of a constant minimum uncertainty 
h are not trivial. A priori, without any hypotheses, the results of the present 
work are meaningful only for classical systems. 

The idea of classical uncertainty relations can be seen in the interpre- 
tation of quantum mechanics on the basis of the hypothesis of Cantorian 
(fractal) space-time [27, 28, 29, 30, 31]. Some authors discussed the quantum 
uncertainty relations due to fractal quantum paths[28, 29, 32, 33]. Their 
reasoning is in fact extendable to macroscopic non quantum fractal paths as 
discussed in [28, 29]. This is the case of classical uncertainty relations arising 
from fractal geometry. 

As a matter of fact, a possible classical background of the Heisenberg's 
relations was proposed long ago by Furth[25] and then discussed by many 
auhtors[24, 26]. There the authors investigated an uncertainty relation of 
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Brownian motion for position and drift velocity of diffusion, the latter is 
not the instantaneous velocity (or momentum), i.e., the conjugate variable 
of position used in Eq.(24). It was believed that [26] the Bohr's concept of 
complementarity also enters in classical mechanics, but with a statistical av- 
erage of momentum; this implies that there is no classical analogs of quantum 
commutation relations. The present work shows that the classical analogs of 
quantum relations exist. 
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Figure 1 




Figure 1: An example of random dynamics: the diffusion of scent. At time 

t a , the molecules of the perfume get out of the bottle at point a. At time 

tf,, the molecules arrive at different points bi (i = 1,2,...) (first dynamical 

uncertainty). For a pony putting his nose at a given point b, all the molecules 

it receives at time tb may arrive there via different paths (second dynamical 

uncertainty). Obviously, as shown by the dotted lines in the figure, the 

second uncertainty cannot take place without the first one. 
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